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Abstract. We explore the relevance of the color-flavor locking phase in the equation of state (EoS) built 
with the Field Correlator Method (FCM) for the description of the quark matter core of hybrid stars. 
For the hadronic phase, we use the microscopic Brueckner-Hartree-Fock (BHF) many-body theory, and 
its relativistic counterpart, i.e. the Dirac-Brueckner (DBHF). We find that the main features of the phase 
transition are directly related to the values of the quark-antiquark potential Vi, the gluon condensate G 2 
and the color-flavor superconducting gap A. We confirm that the mapping between the FCM and the 
CSS (constant speed of sound) parameterization holds true even in the case of paired quark matter. The 
inclusion of hyperons in the hadronic phase and its effect on the mass-radius relation of hybrid stars is also 
investigated. 

PACS. 25.75.Nq Quark deconfinement, quark-gluon plasma production, and phase transitions - 26.60.-c 
Nuclear matter aspects of neutron stars - 97.60.Jd Neutron stars 


1 Introduction 

The possible appearance of quark matter (QM) in the in¬ 
terior of massive neutron stars (NS) is currently one of 
the main theoretical issues in the physics of compact stars 
[T]. The existence of two NS of about two solar masses has 
been confirmed by recent observations m . Based on a 
microscopic nucleonic equation of state (EoS), one expects 
that in such heavy NS the central particle density reaches 
values larger than l//m^, where in fact quark degrees of 
freedom are expected to appear at a macroscopic level. 
Unfortunately, while the microscopic theory of the nucle¬ 
onic EoS has reached a high degree of sophistication [1I51 
[M71I5) . the QM EoS is still poorly known at zero tempera¬ 
ture and at the high baryonic density appropriate for NS. 
In fact the essential theoretical tool, i.e. lattice formulation 
of the quantum chromodynamics (QCD) is inapplicable at 
large baryon densities and small temperature due to the 
so-called Sign Problem [9], and this is due to its compli¬ 
cated nonlinear and nonperturbative nature. On the other 
hand, in the large temperature and small density region 
lattice QCD simulations have provided controlled results 
for the EoS as well as for the nature of the transition [M 

m- 

The mass of a NS can be calculated by solving the 
Tolman-Oppenheimer-Volkoff (TOV) equations with the 
relevant EoS as input. The hybrid EoS including both 
hadronic matter and QM is usually obtained by combin¬ 
ing EoSs of hadronic matter and QM within individual 
theories/models. 


Send offprint requests to-. 


Continuing a set of investigations using different quark 
models [Ullllllllillll], recently we have explored the 
nature of the phase transition with the Field Correlator 
Method (FCM) model of quark matter [TTlIlSlIT^ . which 
in principle is able to cover the full temperature-chemical 
potential plane. In our previous papers we tested 

the FCM model by comparing the results for the neutron 
star masses with the existing phenomenology, which puts 
strong constraints on the parameters of the model, i.e. 
the quark-antiquark potential Vi and the gluon conden¬ 
sate G 2 . 

Recently, we found that the FCM model can be ex¬ 
pressed in the language of the ” Constant Speed of Sound” 
(CSS) parameterization [^1^ . and we showed how its 
parameters can be mapped on to the CSS parameter space. 
We remind that the CSS scheme is a general parameteriza¬ 
tion suitable for expressing experimental constraints in a 
model-independent way, and for classifying different mod¬ 
els of quark matter and establishing connections among 
them. It is applicable to high-density equations of state 
for which: (a) there is a sharp interface between nuclear 
matter and quark matter, (b) the speed of sound in the 
high-density matter is pressure-independent in the range 
between the first-order transition pressure up to the max¬ 
imum central pressure of neutron stars. Given the nuclear 
matter EoS eNM(p )5 the high-density EoS can be expressed 
as 

, , _ reNM(p) P <Pt vans 

^^NM(Ptrons) H" H" ^QM^P Ptrans) P ^ Ptrans 

( 1 ) 
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where the three parameters: the pressure ptrans at the 
transition, the discontinuity in energy density Ae at the 
transition, and the speed of sound cqm characterize com¬ 
pletely the high-density phase. 

In this work we elaborate more on that point by ex¬ 
tending the FCM model in order to include the color 
superconductivity through the color-flavor locking (CFL) 
mechanism, which mimics an explicit dependence of the 
gluon condensate G 2 on the quark chemical potential, as 
initially studied in ref. m- We find that the value of the 
hybrid star maximum mass depends strongly on the FCM 
model parameters, i.e. Vi, G 2 and the gap A. We also 
find that the mapping on to the CSS parameterization 
still holds true when including the CFL phase, though the 
parameter region explored depends on the value of the 
gap A. We also pay special attention to the analysis of 
the hadron-quark phase transition when hyperons are in¬ 
cluded in the hadronic phase. In fact, given the strong soft¬ 
ening of the EoS due to hyperons, it has been often found 
that no quark matter phase transition can take place [13 

[SI- 

This paper is organized as follows. In the next section 
we briefly review the EoS of the hadronic sector, particu¬ 
larly the BHF and DBHF microscopic approaches. In Sect. 
3 we discuss the quark matter EoS, and in Sect. 3.1 we 
illustrate the FCM at finite density, with the inclusion of 
the color-flavor locking effect in Sec. 3.2. Sect. 4 contains 
numerical results, with some details of the EoS for neu¬ 
tron star matter and the hadron-quark phase transition in 
Sect. 4.1. In Sect. 4.2 we discuss the mass-radius-central 
density relation for hybrid stars, and the FCM mapping 
onto the CSS parameterization. Effects due to the inclu¬ 
sion of hyperons are explored in Sect. 4.3. Finally we draw 
our conclusions in Sect.5. 


2 The BHF and DBHF EoS of Nuclear 
Matter 

Empirical properties of infinite nuclear matter can be cal¬ 
culated using many different theoretical approaches. The 
amount of experimental and observational data obtained 
in the last few years, and the intense theoretical efforts 
aimed at their interpretation, call for a firm microscopic 
approach to the modeling of the Equation of State (EoS). 
In a microscopic approach, the only input required is a 
realistic free nucleon-nucleon (NN) interaction with pa¬ 
rameters htted to NN scattering phase shifts in different 
partial wave channels, and to properties of the deuteron. 
In this paper we adopt the non-relativistic Brueckner- 
Bethe-Goldstone (BBG) method [23] and its relativistic 
counterpart, the Dirac-Brueckner- Hartree-Fock (DBHF) 
approximation |25j . 

The Brueckner-Bethe-Goldstone (BBG) theory is based 
on a linked cluster expansion of the energy per nucleon 
of nuclear matter (see Ref. [23], chapter 1 and references 
therein). The basic ingredient in this many-body approach 
is the Brueckner reaction matrix G, which is the solution 
of the Bethe-Goldstone equation 


1 I \kakb)Q{kakb\ ^ 

G[p]U}]=v+ -^ 7 y^G[p;w], (2) 


kakb 


Lu - e{ka) - e{kb) 


where v is the bare NN interaction, p is the nucleon num¬ 
ber density, and ui the starting energy. The single-particle 
energy e(fc) (assuming 


e{k) = e(fc; p) = — + U{k] p), 
2m 


(3) 


and the Pauli operator Q determine the propagation of 
intermediate baryon pairs. The Brueckner-Hartree-Fock 
(BHF) approximation for the single-particle potential U 
using the continuous choice is 

Uik;p)= {kk'\G[p-eik) + eik')]\kk’)a, (4) 

k' Kkp 

where the subscript “a” indicates antisymmetrization of 
the matrix element. Due to the occurrence of U(k) in 
Eq. (13), the above equations constitute a coupled system 
that has to be solved in a self-consistent manner for several 
momenta of the particles involved, at the considered den¬ 
sities. In the BHF approximation the energy per nucleon 
is 

^ = + ^ ^ (fcfc'|G[p;e(A:)-be(fc')]|fcfc')a- (5) 

^ ^ k,k'<kp 

In this scheme, the only input quantity we need is the 
bare NN interaction v in the Bethe-Goldstone equation 
(1). It has been shown that the nuclear EoS can be calcu¬ 
lated with good accuracy in the Brueckner two hole-line 
approximation with the continuous choice for the single¬ 
particle potential, since the results in this scheme are quite 
close to the calculations which include also the three hole¬ 
line contribution 12 ^ 1271123 . The dependence on the NN 
interaction, also within other many-body approaches, has 
been systematically investigated in Ref.[7]- 

However, it is commonly known that non-relativistic 
calculations, based on purely two-body interactions, do 
not reproduce correctly the saturation point of symmet¬ 
ric nuclear matter, and that three-body forces (TBF) are 
needed to correct this deficiency. For that, TBF are re¬ 
duced to a density dependent two-body force by averag¬ 
ing over the generalized coordinates (position, spin and 
isospin) of the third particle, assuming that the probabil¬ 
ity of having two particles at a given distance is reduced 
according to the two-body correlation function [41I61I23 . 

In this work we will illustrate results obtained using 
a phenomenological approach to the TBF, which is based 
on the so-called Urbana model, and consists of an attrac¬ 
tive term due to two-pion exchange with excitation of an 
intermediate A resonance, and a repulsive phenomenolog¬ 
ical central term [331SI]. Within the BHF approach, those 
TBF produce a shift of about -1-1 MeV in energy and —0.1 
fm“^ in density. This adjustment is obtained by tuning the 
two parameters contained in the TBF, and is performed 
to get an optimal saturation point [IE]. 
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Besides a purely phenomenological model, microscopic 
TBF have also been derived and a tentative approach pro¬ 
posed using the same meson-exchange parameters as the 
underlying NN potential. Results have been obtained with 
the Argonne uis, the Bonn B, and the Nijmegen 93 poten¬ 
tials However, at present the theoretical status of 

microscopically derived TBF is still quite rudimentary, Al¬ 
ternatively, latest nuclear matter calculations mMM used 
a new class of chiral inspired TBF [3S1IM1I3Z], showing 
that the considered TBF models are not able to repro¬ 
duce simultaneously the correct saturation point and the 
properties of three- and four-nucleon systems. 

Recently, it has been shown that the role of TBF is 
greatly reduced if the NN potential is based on a realis¬ 
tic constituent quark model |381I39) which can explain at 
the same time few-nucleon systems and nuclear matter, 
including the observational data on Neutron Stars and 
the experimental data on heavy-ions collisions. Moreover 
it has been found that at the highest densities the three- 
hole-line diagrams can give a contribution larger than the 
two-hole-line diagrams, and this can be related to the char¬ 
acteristic nonlocality of the repulsive core [10] as produced 
by the quark-exchange processes. 

In the past years, the BHF approach has been extended 
in order to include the hyperon degrees of freedom miiia, 
which play an important role in the study of neutron star 
matter. In fact, hyperons are expected to appear in beta- 
stable matter already at relatively low densities of about 
twice nuclear saturation density, thus producing a soften¬ 
ing of the EoS with a strong decrease of the maximum 
mass. If this is the case, the existence of heavy NS would 
question the presence of hyperons in their interior, thus 
requiring alternative scenarios. It is therefore of great im¬ 
portance to carry out accurate theoretical calculations of 
hypernuclear matter starting from the available informa¬ 
tion on both nucleon and hyperon interactions. There exist 
several hyperon-nucleon (NY) potentials fitted to scatter¬ 
ing data, i.e. NSC89 [13], NSC97 jH], and ESC08 |35] . 
while the hyperon-hyperon (YY) potentials have presently 
to be considered rather uncertain or unknown, which is ba¬ 
sically due to the lack of appropriate experimental data. 
An alternative description of the hyperon-nucleon system 
has been recently achieved at next-to-leading order in chi¬ 
ral effective field theory gniiiTj. 

In the Brueckner scheme, we have used the phenomeno¬ 
logical NY potentials [3311331133] as fundamental input, and 
found very low maximum masses of hyperon stars, below 
1.4 Mq {Mq = 2 X lO^^g). A proposed solution of this 
so-called hyperon puzzle focuses on the role played by hy- 
peronic three-body forces, and several attempts have been 
made in this direction [331133113(71133] . However, many in¬ 
consistencies still remain, and the solution to this problem 
is still far from being understood. 

The relativistic framework is the one on which the 
nuclear EoS should be ultimately based. The best rela¬ 
tivistic treatment developed so far is the Dirac-Brueckner 
(DBHF) approach [2^. The DBHF method can be de¬ 
veloped in analogy with the non-relativistic case, i.e. the 
nucleon inside the nuclear medium is viewed as a dressed 


particle in consequence of its two-body interaction with 
the surrounding nucleons. The two-body correlations are 
described by introducing the in-medium relativistic G- 
matrix. The DBHF scheme can be formulated as a self- 
consistent problem between the single particle self-energy 
E and the G-matrix. It has been shown that the DBHF 
treatment is equivalent [52] to introducing in the nonrel- 
ativistic BHF the three-body force corresponding to the 
excitation of a nucleon-antinucleon pair, the so-called Z- 
diagram [58] . which is repulsive at all densities, and con¬ 
sequently produces a saturating effect. Actually, including 
in BHF only these particular TBF, one gets results close 
to DBHF calculations, see Ref. [33]. Generally speaking, 
the DBHF gives in general a better saturation point than 
BHF, and the corresponding EoS turns out to be stiffer 
above saturation than the one calculated from the BHF 
-I- TBF method. In the relativistic context the only NN 
potentials which have been developed are the ones of one- 
boson exchange type. In the calculations shown here the 
Bonn A potential is used [3^. Recently, the properties of 
neutron-star matter including hyperons have been investi¬ 
gated within the DBHF approach [33]. In the calculation, 
the effect of negative-energy states of baryons was partly 
taken into account, as well as both time and space compo¬ 
nents of the vector self-energies of baryons and the scalar 
ones. A value of 2.08 Mq was obtained for the maximum 
neutron-star mass, consistent with the recently observed, 
massive neutron stars [3113] . 


3 The quark matter EoS 

3.1 The Field Correlator Method 


The approach we follow to describe the quark matter EoS 
was introduced in [nKunn] ; see ref. [56] for a review. Eor 
our purposes, we are specifically interested in the exten¬ 
sion of this approach to finite baryon density and temper¬ 
ature, and all we need is the expression of the pressure as 
a function of the EoS thermodynamical parameters, i.e. 
the baryon chemical potential /is and the temperature T. 
This is derived in [3711351133] . and below we report its ex¬ 
plicit form . The full pressure, Pqg, is the sum of the gluon, 
Pg, the quark, Pg, and the vacuum, Py, contributions 


Pqg — Pg P^ + Py 


( 6 ) 


j—u,d,s 


where the sum is extended to the three light quark flavors. 
The gluon pressure is 


Q'T'4: poo 


I 


exp(x-b - 1 


(7) 


while the quark pressure for each single flavor with mass 
rUg and chemical potential pLg, is 


rp 4 

A = — 

TT^ 


h^q - Vl/2 


T 


) + 


h-q + 74/2, 


T 


( 8 ) 
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where 


/■°° 1 

= / du -- . -r^, (9) 

Jo Vu^ + (exp — aj + 1) 

being v = mq/T. Finally which represents the pres¬ 
sure difference between the vacua in the deconfined and 
confined phases, is given by 


Pv = 


(11 - iNf) 

32 2 


( 10 ) 


where the number of light flavors in our case is Nf = 3. 

Then, once the quark chemical potentials are related to 
the baryon chemical potential ^b ^ the full pressure Pqg is 
defined in terms of the two parameters Vi and G 2 appear¬ 
ing in Eqs. wm, where Vi indicates the large distance 
static qq potential and G 2 is the gluon condensate. The 
former is essentially of nonperturbative nature and can be 
expressed in terms of an integral of a fundamental QCD 
correlator however there is no direct measure¬ 

ment of its value. The latter is known from QCD sum rules 
PIIHT] . G 2 =< >= 0.012 GeV^ although an 

uncertainty of about 50% affects this estimate. 

It is also interesting to notice that G 2 appears only in 
the vacuum contribution to the pressure, and Py in Eq. 
(fTUl) has the same role of the bag constant of the MIT 
bag model. Moreover, if one turns off the potential Vi, the 
quark pressure Pq becomes the pressure of free quarks, 
and in this case the ECM model reduces to the simplest 
version of the bag model. Therefore Vi can be regarded 
as the main correction to the free quarks dynamics inside 
the bag. 

In addition to the poor knowledge of the phenomeno¬ 
logical values of Vi and G 2 , one has also to deal with the 
dependence of these parameters on the thermodynamical 
variables fis and T. In fact, for the T dependence some 
indications can be obtained from the analysis of the decon¬ 
finement phase transition at T = Ty and fj,B = 0, which 
is supported by lattice calculations. Eor instance, the fact 
that the gluon condensate G 2 is substantially T indepen¬ 
dent except at Tc, where it is sharply reduced by one half 
mm, was already accounted for in the vacuum pressure 
difference of the two phases derived in and reported 
in Eq. (fTUl) . 

As far as Vi is concerning, the following expression re¬ 
lating Vi{Tc), G 2 and Ty is derived within the ECM model, 
in [571155] . 


Tr,= 


aoG^/" 



( 11 ) 


being ag = Then, in |64j it is shown that, 

with G 2 = 0.012 GeV"*^ and with the lattice estimates of 
the critical temperature, Tc = 147±5 MeV or Tc = 154±9 
MeV, Eq. (fTTl) yields Vi(rc) < 0.15 GeV. Another anal¬ 
ysis, [21], based on a fit to the lattice determination of 
the interaction measure (e — 3p)/T^ at several values of 
the temperature around and above Tc, suggests a larger 


potential, Vi{Tc) ^ 0.5 — 0.6 GeV, which, according to 
Eq. (fTTll computed at these Tc, requires a smaller gluon 
condensate, G 2 — 0.003-F 0.004 GeV^. In addition to these 
estimates, an analytic expression for Vi{T), which allows 
to relate Vi (Tc) to the value of the potential at zero tem¬ 
perature, Vi(T = 0), was derived within the ECM in |UT] 
and, in particular, it turns out that Vi(Tc) ~ 0.5 — 0.6 
GeV corresponds to Vi(T = 0) ~ 0.8 4- 0.9 GeV. 

At this point, it is essential to recall that all these re¬ 
sults are obtained at zero baryon density while, in the core 
of NS, densities as large as many times the nuclear mat¬ 
ter saturation density and low temperatures are expected. 
Unfortunately, no lattice simulation can be performed in 
QCD at high ^b, and therefore no numerical indication 
on the density dependence of the two parameters of the 
ECM is available. On the other hand, in the extension of 
the ECM at finite chemical potential discussed in [55] , the 
authors claim that the potential Vi is expected to be in¬ 
dependent oi fiB at least for small values oi ^b- However 
this statement cannot be straightforwardly extended to 
the region of very large density where the environment is 
strongly modified: the number of antiparticles becomes 
much smaller than the number of particles. Therefore, 
even if the interaction strength of quark-antiquark pairs 
is larger than that of quark-quark pairs, the latter plays a 
more important role because of the dominance of quarks 
over antiquarks. In addition, within this framework the 
particle pairing can lead to the appearance of new phases 
that will be discussed in Sect. o 

In view of all the above considerations, it is evident 
that we have very few indications on Vi at large baryon 
density and low temperature, and therefore the best strat¬ 
egy is to treat it as a free parameter of the ECM, with no a 
priori assumed dependence on T or p-B, which effectively 
measures the large distance interaction strength within a 
finite quark density environment, and which has to be con¬ 
strained by the observational data on heaviest NS, through 
our analysis. Accordingly, we do not try to relate Vi(/is) 
to Vi{fiB = 0) and in particular to those values of Vi 
at deconfinement transition, quoted above. Therefore, the 
widest acceptable range of Vi is explored: so, for instance, 
even the case lA < 0 is examined, although, as we shall 
see below, it does not produce sufficiently heavy hybrid 
NS. Finally, as explained more in detail in Sect. 13.21 the 
other parameter of the FGM, G 2 , is also treated a free pa¬ 
rameter in our analysis, essentially for the same reasons 
discussed in the case of Vi. 


3.2 The color-flavor locking effect 

In ref.dl] the dependence of G 2 on qiB was explicitly 
studied by introducing a particular ansatz for G 2 (/iB), 
which was based on previous analysis of the expectation 
of the gluon condensate in dense nuclear matter [5511551 
[57] and in two-color, Nc = 2, quantum chromodynamics 
[551I5U] . In particular, in [551I5U] it is shown that the ex¬ 
plicit form obtained for G 2 (p.b) (it starts as a decreasing 
function at small p,B and turns into an increasing func¬ 
tion at larger piB, with a minimum between these two 
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regions) is related to the effect of diquark pairing and the 
appearance of a mass gap. In fact, at a qualitative level 
this effect presents strong similarities with the occurrence 
of the color-flavor locking (CFL) superconductive mecha¬ 
nism in standard = 3 QCD [7S1IZT], which, as far as 
the total pressure of the quark system is concerned, in¬ 
duces an additional pressure term parameterized by a gap 
A. Clearly, this new pressure term has an essential role in 
shaping the dependence of the gluon condensate G 2 , 
which eventually has the same qualitative features of the 
curve derived in |69) . 

Therefore, due to the unavoidable arbitrariness asso¬ 
ciated to the choice of the ansatz for 6*2 (/rs) encountered 
in m, rather than following that approach we prefer to 
adopt here the same point of view taken for the other FCM 
parameter Vi, i.e. we take G 2 as a free parameter indepen¬ 
dent of T and fiB, that sets the pressure and the energy 
density of the vacuum. Then, the FCM described by the 
two free parameters Vi and G 2 , is suitable to study the 
high density region, and the possible additional contribu¬ 
tion due to the new CFL phase, associated to quark-quark 
pairing, has to be taken into account by adding the addi¬ 
tional CFL pressure contribution to the full FCM pressure 

p 

^qg- 

The presence of color-flavor locked quark matter is ex¬ 
pected at very high pB, and it is realized through quark- 
quark pairing under the constraint that the densities of 
the three flavors, up, down and strange, are equal [nun 
rzHirTi] . The global effect of this pairing on the pressure is 
the presence of the additional term 


Pcfi = 


Stt^ 


( 12 ) 


only when the chemical potential is greater than pb = 
3ms/(4Z\) and the gap A is expected to be in the range 
10 — 100 MeV in the region of interest of pB for the NS. 
Finally the total pressure of the quark matter phase is 
obtained by adding P^fi to Pqg given in Eq. (jO]), and is 
treated as a function of the baryon chemical potential pB 
with three free parameters, namely the potential Vi and 
the gluon condensate ( 72 , coming from the FCM model, 
and the gap Z\, due to the CFL pairing. 


4 Numerical Results 

4.1 EoS of dense matter in beta equilibrium 

In order to study the structure of NS, we have to calcu¬ 
late the composition and the EoS of cold, neutrino-free, 
charge-neutral, and beta-stable matter, characterized by 
two degrees of freedom pB and pe, the baryon and charge 
chemical potentials. The corresponding equations are 

Pi = hpB - qtPe , ^ Piqt = 0 , (13) 

i 

bi and qi denoting baryon number and charge of the par¬ 
ticle species i = n,p,e,p in the hadron phase and i = 
u, d, s, e, p in the quark phase, respectively. 


As far as the hadronic phase is concerning, the Brueck- 
ner calculation yields the energy density of baryon/lepton 
matter as a function of the different partial densities, 

E 

^^Pni Ppi Pei — ^Pn'^n PpUZp) (^Pn Pp)~^i.Pni Pp) 

+ e.fj_{pg) + eel{Pel)i (14) 


where we have used relativistic and ultrarelativistic ap¬ 
proximations for the energy densities of muons and elec¬ 
trons, respectively m- In practice, it is sufficient to com¬ 
pute only the binding energy of symmetric nuclear matter 
and pure neutron matter, since within the BHF approach 
it has been verified [41176] that a parabolic approximation 
for the binding energy of nuclear matter with arbitrary 
proton fraction x = Pp/p, p = Pn + Ppi is well fulfilled, 

^{Pi x)-^{p,x = 0.5) + (1 - 2xfE^y^{p) , (15) 

where the symmetry energy Egym can be expressed in 
terms of the difference of the energy per particle between 
pure neutron (x = 0) and symmetric (cc = 0.5) matter: 


Psyuiip') — 


ld^{E/A) 
8 5a;2 


§(P,0)-^(P)0.5). (16) 


Once the energy density is known (Eq. (fT4ll l. the various 
chemical potentials (of the species i = n^p,e,p) can be 
computed straightforwardly. 


/i* = , (17) 

9Pi 

and the equations for beta-equilibrium o allow one to 
determine the equilibrium composition {pi} at given baryon 
density p and finally the EoS, 


p ^ 2 d e{{p^{p)}) 

^ dp p 


de 

= -e = ppB - e . 
dp 


(18) 


As far as the quark phase is concerning, it is neces¬ 
sary to define the relations among the various pq that 
appear in Eq. (jS]) and the variable pB- For this purpose 
we must distinguish two cases, one with Z\ = 0 and the 
other with A 7 ^ 0. In the first case one has Pcfi = 0 and 
the corresponding EoS for quark matter is determined by 
the conditions of /3-equilibrium and charge neutrality and 
baryon number conservation, as expressed by Eg. (1131) . It 
is then sufficient to express each pq in terms of one single 
variable, namely pB-^n the second case, with A ^ 0, the 
charge neutrality condition is realized in a peculiar way 
m- In fact, CFL pairing occurs if the number densities 
of the three flavors are equal 


Pu — Pd — Ps 


(19) 


which implies vanishing electron density, pe = 0 , in order 
to maintain full charge neutrality. As explained in 
Eqs. (1191) with non vanishing strange quark mass, ms 7 ^ 0, 
are acceptable only if pe 7 ^ 0 , but in any case they allow 
us to determine each single pq in terms oi pb- Therefore, 
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Fig. 1. (Color online)The pressure is displayed as a function of the baryon chemical potential fj.B for the FCM quark matter 
and the pnrely hadronic matter. All calculations for FCM have been performed for G 2 ~ 0.006 GeV"^, and several values of 
Vi have been chosen. The solid curves represent the BHF (green) and DBHF (red) EoS. Each panel shows results for different 
values of the gap A, i.e. 0, 50, and 100 MeV. 




Fig. 2. (Color online)The pressure is displayed vs. the energy density for Vi = 0 (magenta curves) and Vi = 200 (violet) using 
as hadronic EoS the BHE (upper panels) and the DBHE (lower panels). All calculations for FCM have been performed for 
G 2 = 0.006 GeV^. The left, middle and right panels correspond to the values A = 0, 50,100 MeV respectively. 
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each time we consider the case with A 0, we use the 
specific condition in Eqs. (fT^ to express in Eq. ([8]) in 
terms of /is. 

Let us now discuss the main features of the hadron- 
quark phase transition, which we assume to be first-order, 
thus performing the Maxwell construction. Fig. [T] shows 
numerical results for the pressure as a function of the 
baryon chemical potential in the hadronic matter and 
quark matter in beta equilibrium. In particular, the green 
(red) solid curves represent the BHF (DBHF) EoS, whereas 
the remaining curves are the results for the FCM model 
with different choices of the quark-antiquark potential Vi 
(expressed in MeV). For completeness, a negative value 
of the potential, Vi = —50 MeV is also included in this 
analysis. In the left, middle and right panels the value as¬ 
sumed for the gap A is respectively equal to 0, 50 and 
100 MeV. All calculations shown in Fig. [T]are performed 
taking G 2 = 0.006 GeV'^. We notice that with increasing 
the value of Vi the transition point is shifted to larger val¬ 
ues of the chemical potential, hence of the baryon density. 
However, the exact value depends also on the stiffness of 
the hadronic EoS at those densities. In this case, being 
the DBHF EoS stiffer than the BHF, the transition takes 
place at smaller values of the density. We notice that the 
transition point is affected also by the value of the gap 
A, and is shifted toward smaller /is for larger value of 
the gap. We also see that no phase transition occurs for 
negative values of Vi. 

The resulting EoS, for the several cases discussed, is 
displayed in Fig. [51 where one can directly read off the 
phase transition between hadron matter and quark matter 
under the Maxwell construction. We notice that the phase 
transition is allowed only for Ifi > 0, and that the width 
of the plateau is directly related to Z\. In the case Vi = 0 
(magenta curve) the phase transition takes place at very 
low value of the density, the plateau is quite small and the 
pure quark matter phase starts at density about 3 times 
normal nuclear matter density. By increasing Vi the phase 
transition is shifted to larger values of the energy density. 


4.2 Hybrid star structure 


Once the EoS of the hybrid star matter is known, one 
can use the Tolman-Oppenheimer-Volkoff m equations 
for spherically symmetric NS: 


dp 

dr 

dm 

dr 


Gme (I +p/e){l + Airr^p/m) 


1 — 2Gm/r 


= 47rr^ 


( 20 ) 

( 21 ) 


where G is the gravitational constant and m(r) is the en¬ 
closed mass within a radius r. Given a starting energy 
density Cc, one integrates these equations until the sur¬ 
face r = R, and the gravitational mass is obtained by 
Mg = m{R). The EoS needed to solve the TOV equa¬ 
tions is taken from the neutron star matter calculations 
discussed above, and matched with the crust EoS, which 
has been taken from Refs. m 


As is well known, the mass of the NS has a maximum 
value as a function of radius (or central density), above 
which the star is unstable against collapse to a black hole. 
The value of the maximum mass depends on the EoS, so 
that the observation of a mass higher than the maximum 
mass allowed by a given EoS simply rules out that EoS. 
This is illustrated in Fig. |3l where the relation between 
mass and radius (left panel) and central density (right 
panel) in units of the saturation density po = 0.17 fm“^ 
is displayed. Results are plotted for the case in which the 
BHF EoS is used for hadronic matter. In Fig. |3| we plot 
several cases obtained for different values of Vi, G 2 and 
A and in this example the largest value of the maximum 
mass is observed for large values of Vi = 200 MeV, A = 
100 and G 2 = 0.01, and it is compatible with the largest 
mass observed up to now, i.e. (2.01 ± O.O4)M0 in PSR 
J0348-b0432 |3]. 

Recently, it has been shown in ref. m that the FGM 
equation of state can be accurately represented by the 
so-called GSS parameterization. The basic ansatz is that 
a sharp phase transition occurs to a high-density phase, 
where the speed of sound is density-independent. As al¬ 
ready discussed in ref. [55], in all models of nuclear/quark 
matter one can find the four topologies of the mass-radius 
curve for compact stars: the hybrid branch may be con¬ 
nected to the nuclear branch (G), or disconnected (D), or 
both may be present (B) or neither (A). We will discuss 
in detail the FGM mapping onto the GSS parameteriza¬ 
tion in the next subsection. Here we limit ourselves to use 
the same labels in Figs. O and Hjin order to indicate the 
topology of the mass-radius curve, which is strongly re¬ 
lated to the values of the pressure and energy density at 
the transition point, and to the energy density disconti¬ 
nuity. In the framework of the FGM model, the topology 
is related to the chosen values of Vi, G 2 and A. Using 
Fig. [D as a guide, we can obtain various topologies just 
changing Vi, G 2 and A. For example, when combining 
FGM quark matter to the BHF nuclear matter we find 
that, for unpaired quark matter and Vi = —50 MeV, the 
lowest transition point can be obtained when G 2 > 0.006 
GeV"^. In Fig. jSj the corresponding mass-radius relation, 
obtained with G 2 = 0.01 GeW' is displayed by the blue 
dashed line, which exhibits a branch of stable hybrid stars 
disconnected (D) by the hadronic branch. With increas¬ 
ing Vi the transition point moves to larger values of the 
pressure and the energy density, and as a consequence we 
explore regions of the phase diagram where the topology 
changes. For instance, for lA = 0 we can get both (B) 
connected and disconnected hybrid star branches, whereas 
for Vi = 100 MeV connected (G) hybrid star branches are 
present and, for the largest value of Vi = 200 MeV the 
hybrid branch is absent (A). This is clearly shown by a 
cusp in the mass-radius relation, and all configurations 
with radii smaller than the one characterizing the cusp 
are unstable. Therefore only purely nucleonic stars do ex¬ 
ist in this case. However, the stability of those hybrid star 
configurations is related to the modeling of the deconfine¬ 
ment phase transition, as pointed out in ref. [78] where 
the Gibbs construction was used instead of the Maxwell 
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Fig. 3. (Color online) The mass as function of the radius (left panel) and the central density (right panel) is displayed for 
several values of Vi, G 2 and A. The BHF EoS is used for the hadronic phase. The labels (A), (B), (C) and (D) indicate the 
specific topologies of the hybrid star branch. 





Fig. 4. Same as Fig.3, but for the DBHF EoS. 


method. The additional contribution of the CFL pressure 
to the FCM EoS produces only a shift of the transition 
point, and therefore the topology explored can be differ¬ 
ent than the one of the unpaired case, leaving unchanged 
the phase diagram. 

In Fig. 0] we display the mass-radius (left panel) and 
the mass-central density relation (right panel) in the case 
that the EoS used for the hadronic phase is the DBHF. We 
observe a topology similar to the one displayed in Fig. 01 


except the (B) configurations, which do not appear for the 
chosen set of values used for Vi, G 2 and A. 

Finally we comment on the values of the maximum 
mass. In both cases, either BHF or DBHF EoS for the 
hadronic matter, we see that the largest possible values 
of the maximum mass are obtained only for values of 
Vi > 100 MeV, and that only in the DBHF case maximum 
masses well above the observational limit are possible. In 
fact, the heaviest BHF-I-FCM hybrid star has a mass of 
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Fig. 5. The mapping of the FCM quark matter model onto the CSS parameterization. Results are obtained using the BHF 
(upper panels) and DBHF (lower panels) nuclear matter EoS. The green curves are the phase boundaries for the occurrence 
of connected and disconnected hybrid branches. The dashed black line delimit the region yielded by the FCM model. Within 
that region, the symbols give CSS parameter values for FCM quark matter as G 2 is varied at constant Vi (given in MeV). The 
(orange) cross denotes the EoS with the highest ptrana, which gives the heaviest FCM hybrid star. The left, middle, and right 
panels display results obtained with A — 0, 50, and 100 MeV respectively. 


2.03 Mq, and the heaviest DBHF+FCM hybrid star has a 
mass of 2.31 Mq. Those values are indicated by an orange 
cross in Fig.[5l where we display the mapping between the 
FCM and CSS parameters which are, in addition to the 
particular constant value of the speed of sound cqm, the 
ratio of the pressure and energy density in nuclear matter 
at the transition point, ptrans/^trans, and the ratio of the 
energy density discontinuity and the energy density at the 
transition, Aejctrans- In the upper (lower) panels we show 
results for the BHF (DBHF) hadronic EoS, whereas in the 
left, middle and right panels calculations are reported for 
different values of the gap A — 0, 50, 100 MeV respec¬ 
tively. The solid red line shows the threshold value Accrit 
below which there is always a stable hybrid star branch 
connected to the neutron star branch. This critical value 
is given by [TIHIMIIM] 


ACcrit 

^trans 


1 

2 


3 Ptrans 
2 ^trans 


( 22 ) 


and is obtained by performing an expansion in powers of 
the size of the core of the high-density phase. That result is 
analytical and independent on both Cqj^ and the hadronic 
EoS. The solid (dashed) green lines represent the phase 
boundaries for connected and disconnected branches, and 


Table 1. The minimum and maximum values of G 2 (in units 
of GeW^) are shown for different choices of A and Vi. 


A (MeV) 

Vi (MeV) 

/^mtn 

^2 

^max 

0. 

-50. 

0.007 

0.014 


0. 

0.005 

0.012 


100. 

0.003 

0.01 


200. 

0.003 

0.01 

100. 

-50. 

0.009 

0.014 


0. 

0.007 

0.012 


100. 

0.006 

0.01 


200. 

0.003 

0.01 


are obtained for Cqj^ = 1/3(0.28). Those values span the 
range of Cq^ relevant for the FCM, as discussed in [23]. 

The dashed black contour delimit the region accessi¬ 
ble by the FCM calculation. Above that region, the sym¬ 
bols connected by solid lines show the CSS parameteriza¬ 
tion of the FCM quark matter EoS. Along each line we 
keep Vi constant and vary G 2 over the range indicated 
in Table [1] for the two extreme cases Z\ = 0, 100 MeV 
and different values of Vi, for both BHF and DBHF. In 
Fig. [5| Vi varies from -50 MeV up to the maximum value 
at which hybrid star configurations occur, which is indi- 
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Table 2. The total radius R, the radius of the quark core 
Rq, the radius of the hadronic layer Rh and the crust radius 
Rcrust are given for a hybrid star mass M = 2Mq , for different 
choices of the hadronic EoS and A. All radii are given in km. 


EoS 

(MeV) 

R 

Rq 

Rh 

R-crust 

BHF 

0. 

10.37 

0.055 

9.97 

0.345 


100. 

10.44 

0.215 

9.87 

0.355 

DBHF 

0. 

12.78 

1.27 

10.87 

0.640 


100. 

12.72 

2.42 

9.665 

0.635 


cated by an (orange) cross. For the BHF case that value 
is Vi = 240 MeV, G 2 = 0.0024 GeV^ and for the DBHF 
case it is Vi = 255 MeV, G 2 = 0.0019 GeV'^. The vertical 
black dashed lines indicate the parameter regions accessi¬ 
ble by the FGM and consistent with the measurement of 
a M = 2 Mq. Hybrid stars with mass heavier than 2 Mq 
lie on a very small connected branch on the right side of 
the vertical black dashed lines, and cover a small range 
of central pressures, having a very tiny quark core, with 
mass and radius similar to those of the heaviest purely 
hadronic star, as was already discussed in Ref. [53]. For 
completeness, we display in Table [5] the characteristic ra¬ 
dius of a hybrid star with M = 2Mq obtained with BHF 
and DBHF EoS for the hadronic phase and two extreme 
values for Z\ = 0,100 MeV. We chose typical configura¬ 
tions lying on the vertical black lines plotted in Fig.jSJ We 
notice that the radius of the quark core Rq is bigger for 
the stiffest hadronic EoS, being comprised between 1 and 
3 km, whereas for the soft hadronic EoS the quark core 
radius Rq is not larger that a few hundreds meters. In 
both cases the hadronic layer occupies the largest portion 
of the star, and is characterized by a radius Rh of about 
10 km. The crust radius Rcrust is always smaller than 1 
km. 

Moreover we notice in Fig. |S| that along each line of 
constant Vi, ptrans/^trans grows with G 2 , and this can be 
explained by recalling the linear dependence of the quark 
pressure on G 2 in Eq. (nni), so that, at fixed chemical po¬ 
tential, an increase of G 2 lowers the quark pressure, mak¬ 
ing quark matter less favourable, and shifting the transi¬ 
tion point to higher chemical potential or pressure. This 
was already discussed in Ref. [20] for BHF nuclear mat¬ 
ter, and is equally applicable to DBHF nuclear matter. We 
also see in Fig. [S] that the combination of G 2 and Vi moves 
the curves inside the region accessible by FGM which is 
delimited downward by the dashed black line. Fig.[S]shows 
that the introduction of a color-flavor locking effect char¬ 
acterized by a gap A does not change qualitatively the 
gross features of the phase transition, being the topology 
of the hybrid star branch slightly affected. 

4.3 Effects of hyperons on the phase transition 

It is known that the effect of including hyperons in the 
hadronic EoS is to soften the interior of the NS so that it 
becomes difficult to get masses of the stars as heavy as 2 
solar masses. In our analysis, this is confirmed by the pres¬ 
sure obtained for the two parameterizations NSG89 and 


ESG08, which are plotted in Eig. [B] (respectively, dashed 
and dotted black curves) together with the BHF EoS (solid 
black curve) used in Sect. KT\ In fact, the steeper growth 
with /is of the former two curves with respect to the latter 
is an indication of the greater stiffness of the EoS when hy¬ 
perons are neglected. At the same time we notice that the 
NSG89 and ESG08 parameterizations quantitatively give 
very close results. We stress that in this paper we discuss 
hyperon effects obtained only in the BHF approach. 

Then, it is easy to realise that the inclusion of hyper¬ 
ons puts more stringent constraints on the parameters of 
the quark matter EoS, in order to observe a crossing of 
the pressures in the two phases. In Eig. [B]we report four 
examples of quark matter pressure for different choices of 
the parameters (curves (a),(b),(c),(d)). These curves ex¬ 
plicitly show that the parameter Vi is mainly responsible 
for their slope, while A has a much smaller effect and G 2 , 
being an additive constant to the pressure as shown in Eq. 
(O, produces a global upward or downward shift of the 
quark matter pressure. 

Therefore, when going from Vi = — 50 MeV to Vi = 0 
to Vi = 120 MeV, the corresponding curves (a), (b) and 
(c) become less and less steep, and one can observe three 
representative behaviors: (a) shows a phase transition at 
a crossing point below 1100 MeV with any hadronic EoS; 
(b) has the same crossing as in the case (a) and, when 
compared to the BHF EoS, it remains the favoured phase 
at any pB, but when compared to the NSG89 or ESG08 
parameterizations, one observes a second crossing at larger 
Pb] (c), after an interval in which the pressures for the 
hadronic and quark phases are substantially the same, the 
NSG89 or ESG08 curves stay above (c) which, in turn is 
above the BHF curve. 

Therefore, one learns that the quark pressure can ex¬ 
ceed the pressure of hadronic EoS including hyperons only 
at small Vi, typically well below 100 MeV, and one knows 
from the analysis of Sect. l4.2l that smaller Vi correspond to 
NS with smaller masses. To verify this point, we consider 
the cases (a) and (b) (retaining for (b) only the first cross¬ 
ing of the quark and hadronic pressure) and then derive 
the corresponding mass-radius or mass-central density re¬ 
lation which are given in Eig. [7] for the ESG08 EoS (the 
NSG89 case produces almost indistinguishable results). It 
is evident that in the cases (a) and (b), the masses remain 
below 1.5 M 0 . In addition, for (b) we ignored the sec¬ 
ond crossing point above which, in principle, the hadronic 
phase becomes again favourable, but, in any case, this new 
transition would make the NS even softer thus lowering its 
maximum mass and ruling out the possibility of reaching 
2 Mq. 

Before concluding, we reconsider in the other cases the 
procedure followed for (b), namely the derivation of the 
mass-radius relation obtained by systematically consid¬ 
ering only the transition from hadronic to quark matter 
occurring at the lowest value oi pB, while ignoring other 
potential transitions occurring at higher pB- 

By following this procedure, which is somehow justified 
by our poor knowledge of the hyperon interactions, we 
reconsider in detail the various examples of Eig. [B] and. 
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Fig. 6. Pressure vs. baryon chemical potential corresponding to the hadronic EoS’s, including hyperons (dashed and dotted 
black lines) or without hyperons (BHF) (solid black line), and their crossing with the quark matter pressure evaluated with 
four different parameterisations ((a),(b),(c),(d)) lines). 


for convenience, we report in Fig.|8]an enlargement of the 
region at small fis- In fact, from Fig. [5] it is evident that 
the double transition which is observed at a larger scale 
in Fig. |6]for the case (b), now occurs on a smaller scale 
for the cases (c) and (d). These two examples are realized 
with Vi = 120 MeV and, respectively, A — 100 MeV and 
Z\ = 0 and G 2 is tuned to make the quark matter pressure 
almost tangent to the NSC89 or ESC08. While the case 
(c) crosses for the first time the hadronic EoS’s slightly 
above hb ^ 1000 MeV, where the ESC08 and the BHF 
EoS’s are distinguishable, the crossing of curve (d) occurs 
below fiB ^ 980 MeV, before the onset of hyperons. 

In both cases (c) and (d), the second crossing is very 
close, but according to the assumption made we retain 
only the quark matter EoS at larger /is after the first 
crossing. The interesting point is that, as Vi = 120 MeV 
is rather large for these two cases, one expects large max¬ 
imum masses for the corresponding NS. This can be veri¬ 
fied by looking at the right panel of Fig. [71 where after a 
very small plateau, related to the small difference between 
the quark and hadronic pressure along a rather large inter¬ 
val around the transition point, the mass of the NS grows 
above 1.7 Mq for (c) and up to 1.95 Mq for (d) that is rea¬ 
sonably close to the observational constraint of 2 Mq. It 
is remarkable that very similar results are obtained for the 
maximum NS masses in [55] where the NSC89 parameter¬ 


ization is used for hyperons and a sort of QCD corrected 
bag model for quark matter. 


5 Conclusions 

The FCM extension at finite T and fj.B provides us a 
very simple description of the quark dynamics in terms 
of two parameters, namely the gluon condensate G 2 , that 
parametrizes the vacuum pressure and energy density, and 
hence is strictly related to the bag constant of the MIT 
bag model, and the potential Vi , which summarizes the in¬ 
teraction corrections to the free quark and gluon pressure. 
It is then natural to explore the predictions of the FCM 
on the maximum masses (and their associated radii) of NS 
when these parameters are varied. In order to have a more 
complete picture, we include the effect of color supercon¬ 
ductivity through the CFL mechanism, which amounts to 
the addition of a new free energy contribution written in 
terms of the gap A. 

Clearly, any prediction of a quark matter model on 
the structure of a hybrid NS strongly depends on the nu¬ 
clear matter EoS employed and among the large variety 
of nuclear EoS available in the literature, we focused on 
the non-relativistic BHF EoS and its relativistic counter¬ 
part, the DBHF EoS. These are derived within a solid 
microscopic approach, and give us different predictions on 





12 


G. F. Burgio, D. Zappala: Hybrid star structure with the Field Correlator Method 




Fig. 7. Mass-radius (left panel) and mass-central density (right panel) plots of the hybrid stars corresponding to the quark 
matter EoS (a), (b), (c), (d) together with the hyperon parameterization ESC08 of Fig. [6] 


the NS maximum mass. For completeness, we also ana¬ 
lyzed the inclusion of the hyperon degrees of freedom that 
produces a softening of the nuclear matter EoS with the 
consequent reduction of the maximum mass of the NS. 

With this new set of more refined calculations, we con- 
hrm the trend already observed in Enumns], i.e. the max¬ 
imum mass of hybrid stars grows with the two parame¬ 
ters Vi and G 2 while it decreases when A is increased. 
More interestingly, we extend the mapping developed in 
P5] among the parameters of the FCM and those defining 
the CSS parameterization, by displaying the effect of the 
gap A. In fact, from the various panels of Fig. it is evi¬ 
dent that the border of the area of the diagram accessible 
by the FCM (i.e. the dashed black curve) is not sensitive 
to A and it is totally determined in terms of the CSS pa¬ 
rameters. In particular, even the region corresponding to 
configurations associated to hybrid stars with maximum 
mass greater than 2 M©, which is the triangle-like area 
delimited above by the straight solid red line, below by 
the dashed black curve and finally on the left by the al¬ 
most vertical dashed black segment, is only determined in 
terms of the two CSS parameters reported on the x and 
y axes of the diagrams. Therefore one can conclude that 
a particular configuration with mass around or above two 
solar masses can be realized in the FCM by different pairs 
of G 2 and Vi, depending on the specific value assigned 
to A, i.e. the appearance of a color superconducting gap 
can be mimicked by a shift of the other two parameters. 


Therefore, even the mass of the heaviest hybrid star pre¬ 
dicted by the FCM (the orange crosses in Fig. [S]) does not 
correspond to a unique set of G 2 , Vi and A, while, as seen 
in [23], its value strongly depends on the specific choice 
made for the nuclear matter EoS. 


The inclusion of the hyperons induces dramatic changes 
in this picture. In fact, a regular transition from nuclear 
to quark matter with a stable quark phase up to very high 
chemical potential requires a particular tuning of the FCM 
parameters that leads to very low maximum masses, be¬ 
low 1.5 Mq. We have also observed that it is possible to 
find specific sets of the parameters G 2 , Vi and A such that 
the hadronic and quark matter pressure run very close for 
a large range oi ys and, when looking more in detail, one 
observes multiple crossings of these lines, although at large 
/is the phase that includes hyperons is favoured. For com¬ 
pleteness we analyze these cases by retaining only the first 
crossing from the hadronic to the quark matter phase, and 
neglecting the other transitions at higher chemical poten¬ 
tial. In this case it is possible to tune the FCM parameters 
in such a way to raise the maximum mass up to 1.95 Mq, 
although an explanation supporting the strong assump¬ 
tion on the transition is required in order to accept this 
result. 
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